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Abstract. Parkes pulse arrival-time data for 165 radio pulsars spanning from 1990 to 2011 have 
been searched for period glitches. Forty-six events out of the detected 107 glitches were found 
to be new contributions to the entire glitch population which currently contains approximately 
> '. 400 events. 
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1. Introduction 



Two types of timing irregularities have been hampering a full interpretation of the 
rotation of neutron stars: timing noise - a random process seen in phase residuals and 
glitches - a discontinuous increase in the pulse frequency v. Glitches are transient events; 
Q-f the most stringent upper-bound known for the timescale of the rising edge of v, 40 s, 
was obtained by Dodson, McCulloch & Lewis (2002) in a high time-resolution observing 
program for the Vela pulsar. The frequency jump caused by a glitch is actually small with 
the detected maximum relative size of Af g /f ~ 10~ 5 (Yuan et al. 2010; Manchester & 
Hobbs 2011). After a glitch, the slow-down rate |z> often partially recovers exponentially, 
followed by a long-term linear decrease; permanent changes in v and its time derivatives 
t— I ■ are usually left. These widely observed features are used to model glitches. 
>\ 

2. Observations and Results 

o ■ 

■ Approximately 200 Galactic pulsars have been being observed regularly from 1990 up 

until now by the 64-m radio telescope at the Parkes Observatory located in western NSW, 
Australia. Fig. [1] shows 165 pulsars that are in the sample on the P — P diagram. Each 
£SJ , pulse time-of-arrival (TOA) for a particular pulsar was obtained every two to four weeks. 
t— I ■ The data spans for these pulsars range from 5.3 to 20.8yr with a total of ~1911yr. By 
\ carrying out off-line data reduction, 107 glitches were identified in 36 pulsars and 46 
events have not previously been published. Moreover, the recovery fraction Q and the 
timescale Td were measured for 27 exponential decays; post-glitch v were also well mea- 
sured for those pulsars that exhibit long-term linear increase in v between two adjacent 
glitches. More detailed results are being published in Yu et al. (2012). 



3. Discussion 

Young pulsars have been thought to be appropriate samples to study neutron-star 
interiors. It was soon after the discovery of the first glitch that led to the realisation 
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Figure 1. P - P diagram showing the 165 pulsars in the sample. The various sym- 
bols indicate the pulsars where no glitch has been reported (x), glitches were detected 
within our data span (A) and glitches were previously reported before our observations 
(O)- Data for the pulse period and period derivative are from the ATNF Pulsar Catalogue 
I http://www.atnf.csiro.au/research/pulsar/psrcat/ ; version 1.43). 



that neutron stars are two-component rotators: a rigid-body bulk plus a faster-rotating 
neutron superfiuid (Baym et al. 1969). After twenty- year-observation at Parkes, we now 
have further found i) the bimodal distribution of the glitch fractional size Ai/g/v has been 
even clearer (upper panel, Fig. [5]), ii) glitches have frequently occurred in the pulsars that 
have characteristic ages around 10 kyr where, furthermore, large glitches were generally 
seen, iii) exponential recoveries could be resolved into multiple components correspond- 
ing to different timescales and iv) post-glitch v values are often positive, are generally 
larger than the prediction from the magnetic-dipole model and are proportional to the 
ratio between the slow-down rate and the inter-glitch interval with a proportionality con- 
stant ~10 -3 . These observations have further probed that neutron stars may suffer the 
re-configuration of the crustal plate tectonics and/or the substantial release of pinned 
vortices and, the rotational equilibrium after a glitch established via vortex-drifting may 
occur in multiple regimes resulting in the observed exponential and linear recoveries. 

Never-the-less, the available observations are still not sufficient to entirely interpret 
the glitch phenomenon. As shown in the lower panel of Fig. [3J the observed fraction Q 
of the glitch that recovers exponentially has exhibited a bimodal distribution with two 
peaks locating at approximately 0.01 and 1.0 respectively. For the two-component model, 
the conservation of angular momentum indicates Q = I a /I, where / is the total moment 
of inertia of the neutron star and / s the inertial moment of the superfiuid that gradually 
couples to the effective crust following a glitch. The values of Q ~ 1 might have implied 
that, for some glitches, I s involved most of the neutron superfiuid, though it has been 
shown that the core superfiuid may rigidly couple to the realistic crust in a timescale 
< seconds (Alpar, Langer & Sauls 1984; Alpar & Sauls 1988). Moreover, a very large 
Q ~ 8.7 was detected in PSR J1846— 0258 recently (Livingstone et al. 2010) which is an 
even more complicated situation. In addition, peculiar features have also been observed. 
For instance, quasi-periodic oscillation in timing residuals following a glitch has both 
been seen in PSR B2334+61 (Yuan et al. 2010) and the Vela pulsar where, moreover, 
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Figure 2. Histograms for the fractional glitch size Au g /u (upper panel) and the recovery fraction 
Q of exponential decays (lower panel). In each plot, the data for the blank bars are from the 
literature, whereas those for the shaded bars are from Yu et al. (2012). 

excess delays were observed for the pulse arrival times measured at a lower frequency 
(McCulloch et al. 1990). Further exploration on glitch- related phenomena may require 
intensive observations perhaps with high precision and high observing cadence. 
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